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A POSTERIORI ERROR ESTIMATION IN A EINITE ELEMENT 
METHOD FOR RECONSTRUCTION OF DIELECTRIC 
PERMITTIVITY 

JOHN BONDESTAM MALMBERG 


Abstract. We present a posteriori error estimates for finite element approxi¬ 
mations in a minimization approach to a coefficient inverse problem. The prob¬ 
lem is that of reconstructing the dielectric permittivity e = £(x), x E C 
from boundary measurements of the electric field. The electric field is related 
to the permittivity via Maxwell’s equations. The reconstruction procedure is 
based on minimization of a Tikhonov functional where the permittivity, the 
electric field and a Lagrangian multiplier function are approximated by peice- 
wise polynomials. Our main result is an estimate for the difference between the 
computed coefficient and the true minimizer e, in terms of the computed 
functions. 


1. Introduction 

In this note we study an adaptive finite element method for the reconstruction 
of a dielectric permittivity function e = e(x), x G fl, where C is a bounded 
domain with (piecewise) smooth boundary F. This is a coefficient inverse problem 
(CIP) for Maxwell’s equations, where the dielectric permittivity function e, acting 
as the coefficient in the equations, characterizes an inhomogeneous, isotropic, non¬ 
magnetic, non-conductive medium in Possible applications include detection of 
explosives in airport security and detection of land mines. 

The method studied is based on minimization of a Tikhonov functional, where 
the functions involved are approximated by piecewise polynomials. It is intended 
as a second stage in a two-stage numerical procedure for the reconstruction of 
a dielectric permittivity. On the first stage, described in a good initial 

approximation Eq of the dielectric permittivity function is obtained by a globally 
convergent method. This initial approximation is then refined on the second stage. 

The version of the second stage considered here was introduced in [l^ . Another 
version was studied theoretically and numerically in Si 01- There were two 
main reasons for introducing the new version of the second stage in [l9| . The first 
reason was to handle a discrepancy between theory and implementation which was 
present in the previous version. This discrepancy was primarily due to the fact 
that the dielectric permittivity was approximated by a piecewise constant function, 
while the theory required higher regularity. In spite of that discrepancy, reasonable 
reconstructions were obtained, but it remained to be seen whether the new version 
of could produce even more accurate reconstructions. 

The second reason to introduce the version of [13 was to incorporate the di¬ 
vergence free condition for the electric displacement directly into the differential 


Date: February 27, 2015. 


1 



2 


JOHN BONDESTAM MALMBERG 


equation, without having to introduce an additional stabilizing penalty term as 
was done in 

In Q, an a posteriori error estimate for a Lagrangian functional was derived. 
A similar estimate was given in [l9| . but there the amount of detail provided in 
the proof was, for the sake of brevity, kept to a minimum. Here we give the fully 
detailed proof of that estimate. Moreover, we extend the error analysis also to 
include a posteriori error estimation for the Tikhonov functional, as well as for 
the permittivity function itself. The arguments which we use here could easily be 
adapted to obtain such estimates also for the original version of the second stage 
considered in 

The remaining part of this note is structured as follows: In the next section 
we present the mathematical formulations of the direct and inverse problems and 
present the basic results prior to discretization of the problems. In Section[3]we state 
the finite element formulations and perform the error analysis. Some concluding 
remarks are given in Section |4l 


2. The direct and inverse problems 

Before proceeding with the mathematical statement of the problem, we introduce 
some notation. For the bounded domain H C with boundary T, we write 
:= H X (0, T) and Ft := T x (0, T), where T > 0 is a (sufficiently large) fixed 
time. If A C K", n G N, is a domain, we define the norm H-H^ ^ := ||■||//m(x) and 
corresponding inner product (•, ■)^ ^ := (•, where H^{X) is the L 2 -based 

Sobolev space of order m over X, with respect to the usual Lebesgue measure. To 
simplify notation, we will drop the index m whenever it is zero. 

Let := We define the set of admissible dielectric permittivity func¬ 

tions 


(1) := {w G : 1 < v{x) < Vx G ft, u|r = I, Vu|r = 0} 


for some known but not necessarily small upper bound Emax- The set C/® is defined 
to describe a heterogeneous medium in fl, immersed in a constant background with 
permittivity 1 in \ H. 

Under the assumption that e € we consider Maxwell’s equations for an 
isotropic, non-magnetic, non-conductive medium in U: 


( 2 ) 

( 3 ) 

( 4 ) 


dt 

d(sE) 

dt 


-I- V X E = 0 


-VxH=0 


V • (/jH) = V • (eE) = 0 


in rix, 


in Ur, 
in Ur, 


where H = H(x, t) and E = E(x, t), (x, t) G Ur, denote the magnetic and electric 
fields, respectively, and /i > 0 is the constant magnetic permeability. By scaling, 
we may assume that /r = 1. 

To obtain an equation involving only e and E, we combine the curl of and 
derivative of ([21) with respect to t to obtain the second order equation 
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To incorporate (H]) we proceed as in [l^ to expand V x (V x E) 
and use 


V -E = V- 



V • (eE) Ve ■ E 
e e 


-AE + V(V-E) 


where the term V ■ (eE)/e vanishes in view of ([4]). 

Thus, after completing with boundary and initial conditions, we obtain the sys¬ 
tem 


( 5 ) 


5E 

dv 


Ve-E 


= P 


. X 5E, , 

E(., 0) = —(•,0) = 0 


= 0 


in fir, 


on Tr, 


in n. 


where — u ■ V, u denotes the outward unit normal on T, and P G [L 2 (rT)]^ 
is given Neumann data (see Section 4 of @ for details). For well-posedness of 
problems of this class, we refer to [l^ . 

The mathematical statement of the coefficient inverse problem is: 


Problem 1. Given time-resolved boundary observations G G [L 2 (rT)]^ of the elec¬ 
tric field, determine e G such that E = G on Tt- 


The observations G represents either experimental or (partially) simulated data, 
see 0. 

Uniqueness of the solution of coefficient inverse problems of this type is typically 
obtained via the method of Carleman estimates till. Examples where this method is 
applied to inverse problems for Maxwell’s equations can be found in, for example. 


llOll for simultaneous reconstruction of two coefficients, and 


Qtil 


for bi¬ 


isotropic and anisotropic media. However, this technique requires non-vanishing 
initial conditions for the underlying partial differential equation, which is not the 
case here. Thus, currently, uniqueness of the solution for the problem we study is 
not known. For the purpose of this work, we will assume that unicmeness holds. 
This assumption is justified by the numerical results presented in [8|, l9| . 

We introduce the space U"' := {v £ [H^(Ut)]^ : v(-, 0) = 0} for solutions to the 
direct problem, and := {v £ [iJ^(UT)]^ : v(-, T) = 0} for adjoint solutions. 
Both spaces are equipped with the usual norm and inner product on 
Then, by multiplying the first equation in ([S]) by a test function (f> G and 

integration over fix, we obtain, after integration by parts. 


( 6 ) 


0 = - f + (^^(•, n (•> 0). </>(•- 0))n 

+ (VE, - (f, + (^, V . - (^, z. . 

= - (^f ^ f + YE, + (^, V . 0)^^ - (P, <A)r, 

=: ^(e, E, cf), 


where the second equality holds because T) = 0, ^(^ 0) = 0, ^ = P on Fp, 
and Ve = 0 on F^. This leads to the following weak description of the electric field: 


Problem 2. Given e £ , determine E £ such that ^(e, E, </>) = 0 for every 

cf G . 
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zs{t) 


1 



t 


T-6 T-f 


T 


Figure 1 . Schematic illustration of the cut-off function zs appear¬ 
ing in the Tikhonov functional ([7|) . 


Let Eg G denote the solution to Problem!^ for a given e G 17^. We can then 
define the Tikhonov functional F\U^^ 1R+, 



where a > 0 is a regularization parameter and zg = zs{t) G T]) is a cut-off 

function for the data, dropping from a constant level of 1 to a constant level of 
0 within the small interval {T — 6, T — S/2), S <C T, as schematically shown in 
Figure [H The function zg is introduced to ensure data compatibility in the adjoint 
problem arising in the minimization of 0. 

How to choose the regularization parameter a with respect to the level of noise 
in the data is a widely studied topic. Several methods exist, examples are the 
(generalized) discrepancy principle and iterative methods jl|. For the results 
presented here, we regard a as a fixed parameter. 

As remarked before, the initial approximation Sq is obtained using the globally 
convergent method, as described in, for instance This means in particular that 
if Eq is sufficiently close to an ideal solution e*, corresponding to noiseless data G*, 
and if the regularization parameter a is chosen appropriately with respect to the 
level of noise in the data G, then by Theorem 3.1 of Q, the Tikhonov functional F 
is strongly convex in a neighborhood Af C of Eq- If so, then in particular there 
exists a constant c > 0 such that for every Ei, £2 & Af Ci U^, 


c ||ei - £2\\ve < F'^ei] £1 - £2) - F'(£2; ei - £2), 


( 8 ) 


where F'(e; e) denotes the Frechet derivative of F at e, acting on E. 

Throughout the remaining part of this text we will assume that the hypothesis 
of Theorem 3.1 of [^, and hence strong convexity, holds. Then we may seek a 
minimizer e G I/® of F by applying any gradient based method (such as steepest 
descent, quasi-Newton, or conjugate gradient), starting from Eq. 

Such an approach requires that we compute the Frechet derivative of F, which 
is complicated since it involves the implicit dependence of Eg i 1 e. To simplify 
the analysis, in the spirit of optimal control (see for example [2| | for the general 
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theory and some specific examples), we introduce the Lagrangian 

L{u) := F{e, E) + E, A), 

where u = (e, E, A) G [/ := x x c E := E® x x F{e, E) 

was defined in 0, and Sl{e, E, A) was defined in dH). 

We can now minimize F over C/® by minimizing L over U. With the strong 
convexity as above, this would imply that we solve 


Problem 3. Find u € U such that L'{u\ v) = 0 for every v € V. 


Again we use the notation L'(u] v) for the Frechet derivative of L at u, acting 
on V. It can be shown (see Proposition 1 of 0) that 

r// \ / -\ f 

L (w; v) = —iu] e) + —(u; E) + —(u; A), 
where u = (e, E, A) G 17, u = (e, E, A) G V, and 


— (u; e) ■.= a(e- cq, e)a - (W ’ If’ + ((^ ' ^ (f))n^ ’ 

||(u; E) := ((E - G)z|, E),^ - (ef, f + (VA, VE)^^ 

+ (^V£, E)^^=:^(s, A, E), 

— (u; A) = ^(£, E, A). 


In particular, we note that the solution u = (e, E, A) to Problem [3] must satisfy 
^(s, E, A) = 0 for every A G and jF(e, A, E) = 0 for every E G The 

former means that E solves Problem [5] and the latter that A solves the following 
adjoint problem: 


Problem 4. Given e G , determine A G such that ^(e. A, (/>) = 0 for 

every <f) G 


The functional in Problem 0] was defined in ([9|). The problem can be seen as 
a weak analogue of the following system, adjoint to dH): 

- AA - = 0 in GLt, 

fi\ 

^ = -(E-G)4 on Ft, 

r)\ 

A(-,T) = —(., T) = 0 in a 

These observations will be used in the error analysis to be described below. 
But first we shall make some remarks concerning the relation between the Frechet 
derivative of Tikhonov functional and that of the Lagrangian. 

Let Ue = (e. Eg, A^) be the element of U obtained by taking E^ as the solution 
to Problem0]and A^ as the solution to Problem0]for the given e G . Then, under 
assumption of sufficient stability of the weak solutions E^ and A^ with respect to 
e, the observation that 


F{e) = F{e, E^) = F{e, E^) + ^(e, E„ A,) = L{ue), 
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(as Ee, Ae) = 0) leads to 

f)T 

(10) F'{s; .) = —•)■ 

Estimate (|H1) and identity (1101) will play an important role in the error analysis for 
the Tikhonov functional and for the coefficient. 

3. Finite element formulations and error analysis 

In this section we will give finite element formulations for discretizing Problems[2l 
Oandm After that we will turn to the error analysis. We begin by defining finite¬ 
dimensional analogues of the spaces V^, and V, as well as subsets 

corresponding to and U. 

Let Th '■= {K} be a triangulation of 17 and let Ir be a uniform partition of (0, T) 
into subintervals {tk, tk+i], t/c = kr, k = 0, ..., Nr, of length r = T/Nr- With Th 
we associate a mesh-function h = /i(x) such that /i(x) = diam(Ar) for x e AT G Th- 
On these meshes we define0 



= {v 

G : 

v\k G P^iK) \/K G 

Th}, 


Ul 







= {v 

G F'^'"' 

: v\kxi G [P^{K)]^ 

X P^{I) 

VK GThViG Ir}, 

T/adj 

^h 

= {u 

g yadj 

:v\Kxie[P\K)f 

X P\I) 

ViF eThyie Ir}, 

Vh 

= y^ 

X F^- 

X 



Uh 

= ui 

X F.'^- 

X 




where P”(A) denotes the space of polynomials of degree at most n G N over X, 
and the degree q used in the finite-dimensional analogue of F® is at least 1. 

Using these spaces we can state finite element versions of Problems [2] and S] as 
Problem [5] and Problem [6l respectively, as follows: 

Problem 5. Given e G U‘^, determine E/j G F^'“' such that S’{e, (t>h) = 0 for 

every G . 

Problem 6. Given e G , determine \h G such that sTie, \h, 4>h) = 0 for 
every cf G F^*''. 

The finite-dimensional analogue for Problem [3] is: 

Problem 7. Find Uh = {sh, E^, \h) G Uh such that L'{uh, v) = 0 for every 

V &Vh- 

The same remark that was made in conjunction with Problem |3] is also valid 
here: it holds that E/i solves Problem [S] and \h solves Problem [5] for e = £h- 

We will now focus on estimations of the difference between the solution to Prob¬ 
lem [3] and Problem [T] We begin by introducing some additional notation. For 

V = (e, E, A) G F we denote (with some slight abuse of notation) its interpolant in 
Vh by 

WhV = (n?ie, II/iE, II/iA), 

^Observe that the dependence on the step size r in time is not explicitly included in the 
notation for the finite-dimensional spaces. This is justified by the fact that r should be selected 
with regard to h in accordance with the Courant-Friedrichs-Lewy condition. 
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and the interpolation residual by 


rhV = v - XlhV = {rne, r,jE, rhX). 


We will also need to consider jumps of discontinuous functions over Th and Ir- Let 
Ki, K2 &Th such that dKi n 8X2 = e ^ 0. For x e e we define 




lim u(y) + lim v(y), 
y->x,yeKi y->x,yeA;2 


so that in particular if u = wi', where w is piecewise constant on Th and u is the 
outward unit normal, then {w}g = = {wu)\ki + (wu)\k 2 is the normal jump 

across e. We extend {-jg to every edge in Th by defining {u}g (x) = 0 for x e iFnF, 
K € Th- The corresponding maximal jump is defined by 


Hs W := max |{z;}Jy)|, :>^€K€Th. 

y^oK 


For jumps in time, we define 


( 4 ) := 


( lim {v{tk + s) - v(tk - s)), k = I, ..., Nr - I, 


k = 0, Nr, 


and 


Ht W := max{|{i;}t (4)| , (4+i)|} t e (4, 4+i). 

In the theorems and proofs to be presented, we will frequently use the symbols 
ss and < to denote approximate equality and inequality, respectively, where higher 
order terms (with respect to mesh-size or errors) are neglected. 

We are now ready to present the first a posteriori error estimate, an estimate for 
the Lagrangian. The theorem was first presented in 0 , but with only a very brief 
proof. We will here give the full details of the proof. Let us start by recalling the 
theorem: 

Theorem 1 . (A posteriori error estimate for the Lagrangian.) Let u = (e, E, A) G 
U be the solution to Problem[^ and Uh = (sh, E^,, \h) G Uh be the solution to 
Problem^ Then there exists a constant C, which does not depend on u, Uh, h, or 
T, such that 


\Liu)-Liuh)\<c[{\R,\,h\[^]J)^ 


+ {Rxm, +^|[^]sl)aT 
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where 


R, = »(£» - £o) - / () ■ <) 


di 


^A,n = —e^i' 
dXh 


Rx,r = 


du 


V • E;t(-, t)V ■ A;t(-, t) 
eh 

_ L J £ 

T ^ 2h 

(Eh-G)-^ 


d.i -f- 


h£h 


dt, 


[ OAh ] 

L at Jt 


V-A, 


Sh 


-^Sh, 


J) 


^E, a = —Sh 


r aSh ] r aEh ] 

L at Jt , L ai^ J 


T 2h 

[(ysh-Eh)i^]^ 

2heh 


Ven-Eh^ J^.Eh + Jl^Vsh 

- 2 -- 

4. eh 


i?E.r = ?^-P. 
ah' 

Here Jve^ and Je^ denotes the Jacobi matrices ofS/eh and Eh, respectively, and 
(•)^ denotes matrix transpose. 

Note that if eh is piecewise linear, Jvehlif = 0 for every K G Th, hence the 
corresponding term in i?E, n vanishes in that case. 

In the following proof, and thereafter, C is used to denote various constants of 
moderate size which do not depend on u, Uh, h, or t. 

Proof. Using the definition of the Frechet derivative we get 

L{u) - L{uh) = L'{uh; u-uh) + o{\\u - uh\\y) 

The split u — Uh = (u — IlttM) + (11;,^ — Uh) = rhU + {JlhU — Uh) now gives 
L{u) - L{uh) = L'{uh; rhU + {UhU - Uh)) + o{\\u - Uh\\y) 

= L'{uh\ rhu) + L'{uh-, EhU - Uh) + o(||m - Uh\\y). 

The second term vanishes since TlhU — Uh € Vh and Uh solves Problem [7l and we 
neglect the remainder term o(||m — zih||y) since it is of higher order with respect to 
the error. We are then left with 

L{u) - L{uh) w L'{uh-, rhu) = ^{uh\ rhe) + ^{uh] rhE) + ^{uh; rhX), 

oe aE aX 

and individual estimation of these three terms will give the stated result. 

Starting with the first term, we observe that 
dL 

— {uh; rhe) = a(eh- £o, rhe)^ - 


+ 


(V.A,)E,, V(^) 


We aim at lifting all derivatives from the interpolation residuals, thus we split 
the inner product over Ur in the last term above into the sum of inner products 
over Kt := K x (0, T), K G Th'. 

((V-A,.)E„v(=Sf)) = i: ((V.A.)E,,v(!m)) , 

K£Th 
















ERROR ESTIMATION IN RECONSTRUCTION OF DIELECTRIC PERMITTIVITY 


9 


We now integrate by parts, using the notation DKt ■= OK x (0, T), dK'j, := 
{dK \ r) X (0, T), dK^ := (dK n T) x (0, T), K € %: 


E ((V-V)Ek,v(¥) 


KGTh 


Kt 


Ken 


Ken 


/ V(V-Afa)-Eh 
Eh 




Kt 

/ (V-Xh){V-Eh) 


dK-i 


rh£ 




Kt ' ^ Kt 

(V-Afe)(tA-Efa) 


dK' 


Eh 


rh£ 


dK',, 


-E< 


Ken 


+ E< 


Ken 



Kt 


dK' 


(V-Ah)(V-Efe) 


, rts 


£It 


(V-A^)(ia.E;.) 


We observe that V(V-A/j) = 0 on every Kt, K £ Th, since \h is piecewise linear, 
and that e = 1 on F so that r/ie|r = 0. With this in mind, the above calculations 
yields 


((V ^ A,)E., V v) 

Ken 

In order to obtain a residual defined in the whole of fl, as opposed to one con¬ 
taining terms defined only on edges of elements K £ Th, 'vre should manipulate the 
last term in the above expression further. Observe that 

E (i(^^V)(^'E,),r„E) =1 y: (i{(V.AU(>'^E,)).,r..) , 

Ken Ken ^ 

where the factor ^ appears since every internal edge is counted exactly twice in the 
sum over all elements K £ Th- 
Using the approximation 



[ fdSe, [ Xdx 

JdK JK i^K 

where / denotes the maximum of / over dK (see for instance [3), we finally get 

E (it(^'V)(>''E.)).,.Ae) » j;: (j^[(V.A,)(,..E,)].,r,E) 

Ken ^ Ken ^ 
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which gives 
dL 

The) =a{eh- £o, rhs)^ + rhs)^^ 

- + (rir l(V ■ A.)(. ■ E„)|., 

= {Re, rne)^ . 

We can now estimate r^e in terms of £h, using standard interpolation techniques 
(see for instance 0), as 


\rh£\ <Ch^\D^e\ ^ Ch^ 


r deh 1 
L Ol/ j £ 


= Ch\[^U, 


where denotes derivatives of second order with respect to x. Thus 


de 


{uh\ rue) 


< 


C{\Re\,h\[^]J)^. 


We continue with 
f)T 

— (u,; r,E) = ((E, - 


aE' 


+ (VA,, + (^Vs,, t.e)^ 


Again, we seek to lift derivatives from the interpolation residuals, and so we use 
integration by parts to get 


dXh drhE 


/ ^ drhE \ - W / /p. _ _, 

Qt ’ dt /Qt “ Qt , gt /Q 

k=l 


)n 


_ {^h^, ThE) dt 

fe=l 


= -E 

w 


fc=i 
N^-l 


k^l 


r\2 A 

where we have used the fact that = 0 on each subinterval (tfe-i, tfe), for the 
piecewise linear function Xh- 

Since r/tE(-, 0) = T) = 0, this leaves us with 


7 V ,-1 






We now approximate the boundary terms by terms defined on the whole interval, 
using 

1 


/( 4 ) 


fit) dt, 


' — 1 
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that is 


N^-1 


1 ptk 


^h'E)Jt=tk ~ ^ I f {Sh [^]t, rhE}^ dt 

fc=i fe=i 

= <^ [%'],.’-"EW- 

Moving on to 

{yXh, Vr,,E)5^^ = ^ {VXh, Vr^E)^^ , 

Ken 

we integrate by parts and use the fact that AXh = 0 in every K ^Th to obtain 

(VA,, Vr,E)^^ = ^ (_(AA,,r•,E)^^ + (f^,r,E)^J 


Ken 




Ken 

Ken 

We again approximate inner products over dK^ by inner products over Kt, so 
that 

(VA., Vr.E)^^ « [^]^, r,E)^^ + (^, r,E)^^ . 

Together with previous calculations, this gives 
f)T 

—Xu,; r,E) ((E, - G)^,^ r,E)^^ - [^]^, r,E)^^ 

= {Rxm, + (i?A,r, ■riiE)p^ . 

We once more use interpolation estimates 

la^Ei 


9E' 


|r;,E| < C £i^E + 


dt^ 


)-C(),|[lM.|+r|[it-].|). 


to get 


dL, 

r,-E) 


< 


+ <l«x,rl,k|[fM,|+Hm.l> 


Tt 


It remains to estimate 
BT 

— {un; r,X) = - + (VE,, Vr,X)^^ 

+ Vr,A)^^-(P, r,A)p^. 

Just as before, we obtain 

/g. drhX \ ~ /Ef. ra^l p \\ 

nh Qt , Qt ~ \ T L at Jt> 

and 

(VE,, Vr,A)^^ « [^]^, r,X}^^ + r,X)^^ . 
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Consider the term 

Integration by parts yields 


(It ^ \ 
K€Th 


Kt 


vt„a; 


= (- (V (^). 


RGTh 


= E 


'^Veu + 


/ Kt 




Ken 

/vEh-EhVT_ + Vsh 


= - 


/ SJt 




ReT;. 


where for the second equality we have used the identity V(Ve;i ■ E/i) = + 

^Eh'^^h- 

Noting that V£ft,|r = 0 as £?i G U^, and using the usual approximation for {-jg 
inside elements K G Th get 


(^.VTa)„^»(^Va 




, ThX 


+ {^i(y^h-'EhM,, rh\)^ 


Combining the results for ^{uh', J’/iA) and estimating r?iA in terms of \h just 
as ThE was estimated in terms of Eh gives 


rh\) 


Qt 


<C({|flE,nl.M[w].l+H[T].l>, 

<l«.=.ri.M[fM.|+q[l^].|)J. 


+ 


which completes the proof. 


□ 


One should note that the terms in the error estimate of Theorem [T] which are 
derived from ^{uh', rhX) and §^{uh; ruE) estimate how accurately the solutions 
of Problem [5] and Problem 0] are approximated by the solutions of Problem [S] and 
Problem |6l respectively, for the approximate coefficient £h ■ The remaining term, 
{Re, h I can be interpreted as the error induced by approximating e by 

£h- Thus, if we are mainly interested in that error, or if we can postulate that 
the finite element approximations Eh and A^ are computed with relatively high 
accuracy, then |i?£| may be used as an error indicator by itself. The significance of 
Re will be further illustrated by the error estimates for the coefficient and for the 
Tikhonov functional. 

We now proceed with an error estimate for the coefficient itself. An error estimate 
for the Tikhonov functional will follow as a corollary. 
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Theorem 2. fA posteriori error estimate for the coefficient.) Suppose that the ini¬ 
tial approximation eq and the regularization parameter a are such that the strong 
convexity estimate ([5]) holds. Let u = (e, E, A) S t7 be the solution to Problem\^ 
and let uu = [ek, E/i, \h) G Uh be the solution to Problem^ computed on meshes 
Th and X^. Denote by E and A the solutions to Problem and Problem^ re¬ 
spectively, with permittivity Eh, and set u = (Eh, E, A) G U. Then there exists a 
constant C, which does not depend on u, Uh, h, or t, such that 


Ik ^hWye ^ C(cer]-\- \\R,;\\q), 

where c^ '■= max{l, ||Vc;j|k (q)} and p = p(uh) is defined by 

Proof. Using strong convexity ([8]), we obtain 

Ik - £h\\ve < c{F'{e; E-Eh) - F'ieh; £ - £h)) ■ 
Since 6 minimizes F{e) we have F'{e; e — eh) = 0 and thus 
|2 


Ht 


(11) ||e - e/illye < c|F'(e/i; £ - £/i)| = c 

where the last equality follows from m- 
We expand 


BL 


dL 


dL dL dL 

— (u; E - Eh) - -^{Uh] E-Eh) + -^(Uh; £ - £h) 
OE OE OE 


( 12 ) 


< 


dL \ \ 

— (u; E-Eh) - -^{Uh; £ - £h) 
OS oe 


dL 

— (Uft,; £ - Eh) 


=■■ |0l| + |02|, 


and estimate the two terms |0i| and | 02 | separately. 
For 01 we use the linearization 


01 = 


dmE 

d'^L 


+ 


dXdE 


(uh', E — E/i; £ — c/j) + o(||E — E/i||n.j,^ i) 
(uh', A — \h',£ — £h) + o(||A — A^lloy^ i), 


o2 T T 

where and denote mixed second partial Frechet derivatives of L. Again, 
the remainder terms are neglected as they are of higher order with respect to the 
error. Thus, after exchanging the order of differentiation, we are left with 


(13) 


01 


DAe-e 


dL, ~ ^ , dL, ~ ^ , 
aE E E/i) + A \h) 


where denotes differentiation with respect to the first component in Uh and 

action on £ — £;j. 

We split E — E?i = (E — n?iE) + (IIiiE — E/j) = r/iE + (fl/jE — E?i) and use the 
fact that \h solves Problem S] with coefficient Eh, so that f§(u/t; n?iE — E/j) = 0 
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as n?iE — E/i S This gives 

dL ' -r-, \ 9L ~ dL ^ ~ ^ . dL ~ 

(14) ®“ ^h) = ’^h^) + ^(^^’ 

Similarly, we have 

dL ~ dL ~ dL ~ dL ~ 

(15) ^(^'*’ X-Xh) = -^K^h] rhX) + -^{uh; n/,A - Xh) = rhX). 

as E/j solves Problem [5] with coefficient Sh- 
Combining (IT^ . (fHl) . and m gives 

(£ - + (V (^) r.E, V • A,)^ 




dEh drhX \ 

/nr 


(^) 


Bh, V • r^X 


In the same manner as in the proof of Theorem [1] we integrate by parts in time 
and approximate jumps to get 




_-] Jtk-i 


dt 


fc=i 


{e-eh)rh'E, i 


and 


^ k—1 —1 

1 r^Efc- 


{e-Sh)l [^]t’ 


Thus 


e,< |r-a|i|[l^] 


r/jE 


+ (|v(^)||V-A.|, 


+ (|e - £hl 7 I [^]tl> \ 


J-ftE 


V • ThX 


< Ik - £h\ 


Aoo(n) 


1 Ir OAh ] 
r I L dt J t 




(^)||E‘i’ 


/ s-eh j 

\ ) 


Aoo(f2) 


|V-A;,|, 


r/jE 


IE, 


V • r,A 


Note that 


V ( 

— 

V(d 

-Eh) (E-efe)Veh 

\ eh J 

Lcxd (^) 

e. 


< 

J_ 

Sh 



Aoo(f2) 


+ 


-i^OO (^) 


l|Vs,|| 


ioo(n) 


Ik - ^h\ 


Loo (^) 
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and observe following facts: 




,(n) 


< C ||e - eh\ 


ye ) 


Loo(n) 


< 1, 


P > 0, 


the first following from the Sobolev inequality and the second from noting that 

1 < Cft(x) < Cmax, X e O, by (P). 

Using these observations, and interpolation estimates 


we get 


(16) 




rh'E 


+ r 

s 1 1 

raEfc] 

1 at It 





<c{h\[^] 

+ r 

s 1 1 

raAfc] 

1 at It 

), 


V • 




)> 

c(q 

raxfei 1 
1 at Itl 



11 ^ 


+ l|V£i.|L^(n)(|V-A,|,h 

\dEh] 

L J s 

+ ^l[ 

|\ 

dt J1 1 / Qt 

+ (^ 

1 

11 at It 



-1 1) 
-Itl/ 


+ IIV£.IL^(o)(|e.|, \[^ 

]sl 

I rav-Afc 

II at 

]t )nT) 11^ 

CCeP II 







where Cg and r] were defined in the statement of the theorem. 

Turning to 02 of (HU, we use the techniques of the proof of Theorem [T] to 
estimate 


(17) |02| < C (|i?,|, |c - < C lie - e,||^ < C ||i?,||^ ||e - e„||^. . 

Combining estimates m and dm with (HP and (HP, we conclude that 

l|e-e?i|lve <C{cer]\\e-eh\\ve + ||i?e|ln lk-^/«llvO, 
and the result follows. □ 


Again, just as for the error estimate for the Lagrangian, we see that if the 
numerical errors for solving the direct and adjoint problems are relatively small, 
that is, when u ~ Uh with relatively high accuracy, then ||i?£||Q dominates the error 
estimate. 

Corollary 1. (A posteriori error estimate for the Tikhonov functional.) Under the 
hypothesis of Theorem\^ we have 

\F{e)-F{en)\<c{clr^^ + \m\l), 

with Ce and p as defined in Theorem\M 

Proof. Using the definition of the Frechet derivative and m we get 
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Neglecting the remainder term as it is of higher order with respect to the error, and 
estimating £ — £h) as in the proof of Theorem [21 we obtain 

|J^(e) - F{eh)\ < CiceV + ||i?e||ye) ||£ - ehWye ■ 

Applying Theorem |2] to estimate ||e — e;i||ye, we arrive at 

- F{sh)\ < C{c,r^ + WReWysf < C (cl^f + WR^fy.) . 

□ 


4. Conclusion 

We have presented three a posteriori error estimates for an adaptive finite element 
method for the coefficient inverse problem, Problem [T] for the Lagrangian, for the 
Tikhonov functional and for the coefficient. The latter two are presented here for 
the first time. Each estimator consists essentially of three parts, an estimate for 
the error resulting from finite element approximation of the solution to the direct 
problem, a similar estimate for the finite element approximation of the adjoint 
problem and an estimate corresponding to the approximation of the coefficient. 
The latter part is characterized by the residual in all estimates. 

Explicit solution schemes and numerical testing, including the proper choice of 
regularization parameter a, will be the subject of forthcoming papers. 
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